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^ ■ Multiply gravitationally lensed objects with known time delays can lead to 

direct determinations of Ho independent of the distance ladder if the mass dis- 
tribution of the lens is known. Currently, the double QSO 0957+561 is the only 
lensed object with a precisely known time delay. The largest remaining source 
of systematic error in the H determination results from uncertainty in the mass 
& ■ distribution of the lens which is comprised of a massive galaxy (Gl) and the 

cluster in which it resides. 

We have obtained V-band CCD images from CFHT in order to measure the 
mass distribution in the cluster from its gravitional distorting effect on the ap- 
pearance of background galaxes. We use this data to constuct a two-dimensional 



mass map of the field. A mass peak is detected at the 4.5<r level, offset from, but 
consistent with, the position of Gl. Simple tests reveal no significant substruc- 
ture and the mass distribution is consistent with a spherical cluster. The peak 
in the number density map of bright galaxies is offset from Gl similarly to the 
mass peak. 

We constructed an azimuthally averaged mass profile centered on Gl out to 
2' (AOOh^ 1 kpc). It is consistent with an isothermal mass distribution with a 
small core (r c w 5" = 17/tT 1 kpc). The inferred mass within 1 Mpc is consistent 
with the dynamical mass estimate but 2a higher than the upper limits from a 
ROSAT X-ray study. 

We discuss implications for H in a future paper. 
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1. Introduction 

One of the most interesting cosmological applications of gravitational lensing 
is the determination of Hubble's constant. This can be achieved by measuring the 
time delay between the observed light curves of multiply imaged QSOs (Refsdal, 
1964). If the redshifts of the lens and source are known, and the mass distribution 
of the lens is well-understood, Ho can be determined directly. Two advantages 
of this technique are that it is independent of the distance ladder-based methods 
and it is insensitive to motions of the nearby Universe. 

The first-discovered doubly imaged quasar Q0957+561 (Walsh et al. 1979) 
has been the subject of intensive monitoring campaigns and is currently the only 
lensed system with a well-determined time delay (Thomson & Schild 1994, Schild 
& Thomson 1995, Kundic et al. 1996). The lens is comprised of a primary lens 
(the galaxy Gl) and the cluster containing Gl (Young et al. 1981, Bernstein 
et al. 1993, Grogin & Narayan 1996). The largest remaining uncertainty in Ho 
derived from this system is due to the uncertainty in the mass distribution in 
the lens, specifically, the ratio of the galaxy surface mass density to the cluster 
surface mass density. 

An attempt has been made to measure the mass of Gl spectroscopically, yield- 
ing a velocity dispersion of 300 ±50 km s _1 (Rhee 1991). However, this is a diffi- 
cult measurement due to the proximity of the QSO B image. The spectroscopic 
value is somewhat lower than the ~ 400 km s _1 implied by the Faber- Jackson 
law (Bernstein et al. 1993). 

In this paper we describe the direct measurement of the 0957+561 cluster 
surface mass density by studying its distorting effect on the appearance of back- 
ground galaxies. This technique has recently been demonstrated to work well on 
massive clusters (e.g. Tyson et al. 1990, Tyson & Fischer 1995, Squires et al 
1996), however, the cluster in the 0957+561 field is believed to be of significantly 
lower mass than previously studied clusters and will be a difficult test of the 
method. 

In §0 and §|3| we describe the observations and analysis. §(| describes the 
mass reconstruction techniques. §[5] describes simulations which were carried out 
in order to calibrate the mass measurements and quantify the uncertainties and 
systematic errors in the analysis. discusses the 2-d mass maps and azimuthally 
averaged radial mass profile, §|7] compares the lensing results to previous mass 
determinations and §|| contains the conclusions and a description of some future 
work. 
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2. Observations 

Twenty-nine 720s V-band images were taken of the region around the double 
quasar Q0957+561 using the Canada France Hawaii Telescope on 10-11 Jan 1994. 
FOCAM was used with the thick LORAL3 2048 2 CCD. This chip has 0.207" 
pixels and a field size 6' on a side. The effective field size is reduced somewhat 
because of vignetting at the corners due to the filter wheel. The conditions were 
photometric. The telescope was dithered between observations which allowed us 
to construct a flat-field from our observations. The final combined V-band image 
has a useable field size of 5.8' x 5.5' with FWHM = 0.6" and the RMS noise of 
the sky is approximately 28.7 V mag per square arcsec (Fig. [l]). 

The photometric zeropoint for the 0957+561 V image was determined to an 
estimated accuracy of 0.03 mag from observations of five standard stars in the 
NGC 2264 field (Christian et al. 1985). Reddening in the 0957+561 field is low, 
E(B-V) < 0.03 mag (Burstein & Heiles 1982) so we assume E(B-V) = in this 
paper. 

We located 21 bright isolated stars with V < 23 to investigate the behaviour 
of the point- spread-function (PSF). These had measured ellipticities spanning 
0.00 < e < 0.03 (e = 1-b/a) with a mean of e = 0.015. The 20 objects with 
e > have position angles ranging from —55 <PA (deg) < 10, implying a 
very slight position independent trailing, possibly arising from guiding errors or 
during the image combining procedure. This effect must be very small since any 
uncertainty in the ellipticity measurements of round objects will tend to result 
in an overestimate in the ellipticity on average. We test systematics arising from 
this PSF anisotropy by carrying out simulations using the PSF derived from the 
image (see §Bj). 

3. Faint Galaxy Photometry and Analysis 

The faint galaxy analysis was carried out using the analysis software ProFit 
(developed by PF). This software, starting with the brightest objects, fits an an- 
alytical model to each galaxy, using weighted, non-linear least squares, and sub- 
tracts the galaxy light from the image. It then proceeds to successively fainter ob- 

Fig. 1. — V-band CFHT image of the field centered at 0957+561. The image is 5.8' x 5.5', 
north is at the top and east is to the left. The total exposure time is 5.8 hours. The seeing 
is FWHM = 0.6" and the RMS noise of the sky is approximately 28.7 V mag per square 
arcsec. The image intensity scaling is logarithmic. 
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jects. Once it has detected and subtracted all the objects in an image it replaces 
each galaxy in turn and refits and resubtracts until convergence is achieved. The 
software outputs brightness, orientation, ellipticity and other image parameters 
based on the fitted function. Fig ^ shows the final subtracted image. 

Fig. H shows galaxy counts for the faint galaxies in the field of Q0957+561. 
The magnitudes are isophotal magnitudes with outer isophote of 27.8 mag square 
arcsecond (31.2 mag per pixel). Also shown are recent V-band blank field counts 
from Smail et al. (1995); while the slopes of the counts are in good agreement, 
the Smail et al. counts lie about 45% below or 0.5 magnitudes to the right of our 
counts. Fig. H shows the size-magnitude relationship for the faint galaxies (r^ 
is the radius which contains half the light, uncorrected for seeing). The galaxy 
sequence is well-separated from the stellar sequence down to about 24th mag. 



4. Mass Reconstruction 

For gravitational lensing, the relationship between the tangential shear, 7^, 
and surface mass density, E, is (Miralda-Escude 1991, 1995), 

7r (r) = «(< r) -7e(r), (1) 

where k = H/H crit and r is the radial distance from a given point in the mass 
distribution. The first term on the right is the mean density interior to r and 
the second term is the mean density at r. The basic technique for 2-d surface 
mass density reconstruction using the distortions of faint background galaxies is 
outlined in Kaiser & Squires (1993) (the KS algorithm). For a given coordinate 
on the image (R), the distortion quantity for the i th galaxy is: 



r>f3\ - l-(VoQ 2 [cos(2^)(Ax, 2 - Ayj 2 ) + 2 sin(2^) A^Agj 
A ' ~ 1 + {k/atf X Axf + Ay? ' {Z) 

where {bi/ai) and Q{ are the galaxy axis ratio and position angle, respectively. 
Ax and Ay are the angular horizontal and vertical distances from R to galaxy i, 
and D is related to the tangential shear by (Schneider & Seitz 1995): 

< D(r) >= 2 ^f-f» - (3) 
[1 - K(r)\ 2 + 7j.(r) 



Fig. 2. — The same as Fig. [1] after ProFit galaxy fitting and subtraction. 
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Fig. 3. — Galaxy plus star counts for the observed image (points) and one of the simulations 
(dotted line). The galaxy magnitudes are isophotal magnitudes with outer isophote of 27.8 
mag square arcsecond. The error bars are based on Poisson statistics. The solid line shows 
the input function for stars and galaxies. Also shown are recent V-band counts (dashed line) 
from Smail et al. (1995) which are about 45% below or half a magnitude to right of the 
0957+561 counts. 



- 6- 




Fig. 4. — Half light radius vs. apparent V magnitude for the data (lower panel) and one of 
the Monte-Carlo simulations (upper panel). The radii are not corrected for seeing. Good 
quantitative agreement is seen. The vertical lines indicate the magnitude range of galaxies 
used in the mass reconstruction 
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In the weak lensing regime k « 1, and 7^ << 1, 7^ ~< -D > and the surface 
mass density is given by: 

^=ht^0k?' (4) 

where N is the number of galaxies and n is the number density of galaxies. 
Eqn. |] assumes that the galaxies are intrinsically (in the absence of lensing) 
randomly aligned. The critical density, T, crit depends on the redshift distribution 
of the background galaxies. Because of the random intrinsic alignments of the 
background galaxies, the formal error of k from Eqn. [| is infinite and in practice 
one must employ a smoothing kernel: 

K J nir^i Ax, 2 + Ayi 2 v ; 

In this paper we have have used a smoothing kernel of the form (Seitz & 
Schneider 1995): 

W{ Ax, Ay, s) = 1 - (l + Ax2 +/ y2 ) e"^ 2+A ^ 2 , (6) 

where V is referred to as the smoothing scale. 

The 2-d mass map for the 0957+561 field is shown in Fig. 5. We discuss this 
further in § p.l| . 

Because of the smoothing kernel, plus biases introduced by edge effects in the 
images, this formula is mainly useful for determining the 2-d shapes of mass dis- 
tributions. A less biased way of obtaining mass estimates as well as azimuthally 
averaged density profiles is (Fahlman et al. 1994): 



n _/ \ 1 DAR) 
k[t < n) - K{ri <r <r ) = YTt\ 1^ 



2vrm(l - r 2 Jr 2 ) r ^ ro Axf + Ayf 



(7) 



Fig. 5. — Mass map of the 0957+561 field using the technique of Kaiser & Squires (1993) 
superimposed on a V-band CCD image. A total of 1651 galaxies (arclets) with 24.0 < V < 
26.5 were used in the reconstruction and the the smoothing scale is s = 30". The contours 
are stepped in units of 0.5cx. North is up and East is to the left. The displayed field is 327" 
on a side. 
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r 2 v Di(R) 
2N lo r ^ ro Ax^ + Ay,*' 

where Ni Q is the number of galaxies between r, and ro, and m is the corresponding 
number density. 

In the upper panel of Fig. |6] we show the calibrated 0957+561 cluster radial 
mass density profile resulting from Eqn. [7[ The error bars are based on the la 
galaxy-to-galaxy scatter of the values inside the sum in Eqn. |7|. Details of the 
mass calibration and systematic errors can be found in §[| and a discussion of the 



radial mass profile is in § p72 



It should be mentioned that galaxy distortion is insensitive to flat sheets of 
mass. Consequently, all mass measurements described in this paper are uncertain 
by an unknown additive constant. If there is a substantial flat component to the 
mass distribution our mass estimates will be lower limits (see §6.21). 



5. Calibration and Systematics 

Application of the mass reconstruction techniques described above will result 
in underestimates of the mass due to: measurement error, the effects of seeing, 
and the intrinsic ellipticity of the background galaxies. These effects must be 
calibrated in order to correct the mass estimates. Calibration was carried out 
using two methods. The first involved applying a known constant shear to the 
Hubble Deep Field (HDF) (Williams et al. 1995) and the second was a full 
Monte-Carlo simulation of the data. 



5.1. Hubble Deep Field Simulations 

Compared to the CFHT data the HDF data is both much higher resolution 
and much deeper. Therefore, it offers an excellent opportunity to calibrate our 
mass estimates as described in Kaiser et al. (1995). The technique involves 
stretching the HDF data by 1 + 5, convolving with the PSF and adding noise. 
The values of Di are measured for each galaxy and compared to the unstretched 
values of Di. The quantity of interest is the recovery factor, C =< ADi > J 8. 
The raw values of k are then multiplied by 1/ < C > to yield the corrected 
values. One problem is that the HDF V-band filter (F606W) is significantly 
different from our V-band filter in both central wavelength and width. In order 
to insure that we were sampling the same population of galaxies we chose a 2.5 
magnitude range in the HDF data which yielded the same surface number density 
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Fig. 6. — The upper panel is the radial mass density profile for the 0957+561 cluster from 
Eqn. [7| (weak lensing approximation) centered on Gl. The points are the data for the cluster 
derived from 1307 galaxies having 24.0 < V < 26.5. The solid line is the mean profile and the 
la upper and lower limits for 102 simulations of a spherical cluster having r c = 5" = 17/iT 1 
kpc and S = 36OO/iM pc -2 . r max = 168" = 550/i -1 kpc. The input model is shown by 
the dotted line. The lower panel is the radial density profile from Eqn. || Both the data 
and the simulations have been corrected using the values of k and 7r from the input profile 
for the simulations (see §6.2j). 
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of galaxies as in our field for the range 24.0 < V < 26.5 (see §6.1| ). This turned 
out to be 23.7 < F606W HDF < 26.2. This V-F606W magnitude difference is 
roughly equivalent to the range of values calculated by Fukugita et al. (1995) for 
late-type galaxies at intermediate to high redshift. We stretched the HDF data 
in two different directions and obtained a value of < C >= 0.37 ± 0.03. 

5.2. Monte-Carlo Simulations 

The HDF simulations yield < C > which is appropriate for a constant dis- 
tortion of all galaxies. In reality galaxies at different redshifts are distorted by 
different amounts (i.e. high- z galaxies are more distorted). Therefore, if galaxy 
size in a fixed apparent magnitude range is a function of redshift, one would 
expect C to vary as a function of redshift resulting in a mean value which dif- 
fers from that of the HDF simulations. By carrying out full 3-d Monte-Carlo 
simulations of the observed data we can quantify this difference. Aside from cal- 
ibrating the lensing results, simulations of the observational data are useful for 
quantifying systematics in the inferred mass distributions and determining S cr ^. 

The simulations consist of artificial galaxies distributed in seven redshift shells 
(z = 0.0 - 0.355, 0.355 - 0.5, 0.5 - 0.65, 0.65 - 0.8, 0.8 - 1., 1.0 - 2.0, 2.0 - 7.2). 
For each shell, simulated images are made with pixels equal to one quarter the 
area of the pixels on the observed image. Galaxy images are generated for each 
shell based on the quiescent models of McLeod & Rieke (1995). Briefly, the 
McLeod & Rieke models involve six galaxy types (E, SO, Sab, Sc, Sd, Im) with a 
formation epoch of z = 7.2. H = 60 km s -1 Mpc" 1 and q = 0.05 were assumed. 
Evolution and k-corrections are based on the synthesis models of Bruzual & 
Chariot (1993). Bulge-to-disk ratios are taken from King & Ellis (1985). The 
model counts were normalized to fit our V-band galaxy counts. Stars were also 
added to the simulated images based on the observed star counts (measured for 
V < 24 and extrapolated with the same slope to the detection limit). 

Initially we used galaxy size - absolute magnitude relationships from Im et al. 
(1995), however, these were seen to produce apparent magnitude - size relation- 
ships substantially larger than observed for V > 24. Smail et al. (1995) also saw 
small galaxy sizes for V > 24 with similar image quality. The evolution of galaxy 
sizes (as well as other aspects of galaxy evolution) is not well understood. One 
explanation for the small galaxies we observe is that we are seeing more small 
galaxies at low to intermediate redshift than predicted by the McLeod & Rieke 
models. This might also explain why these models underpredict the density of 
apparently faint galaxies. Alternatively, higher redshift galaxies may have their 
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emission more concentrated towards their centers than local galaxies or we may 
be seeing small emitting regions within the galaxies. It is not currently possible 
to distinguish between these (and other) scenarios. For the present work, we 
arbitrarily assume that the galaxy sizes evolve as (l+z) - ™. We determine n by 
stretching the simulated data by a constant amount as was done for the HDF 
simulations and modifying n until we get the same value of < C > as is observed 
for the HDF data in the range 24.0 < V < 26.5. This yielded n = 0.75; the size- 
apparent magnitude relation for the simulated galaxies are in good agreement 
with the faint galaxy measurements down to the detection limit (see Fig. ^). 

The galaxies in each shell were distorted with various spherically symmet- 
ric mass distributions located at z=0.355 under the assumption that all galaxies 
within a shell lie at the number weighted mean redshift of that shell. The dis- 
torted images were combined and convolved with the slightly elongated PSF 
derived from the observations. The images were block-summed 2x2 and Gaus- 
sian noise was added. The simulated images were processed identically to the 
observed images. 

Fig. |3] shows the observed galaxy counts of one of the simulated data frames 
superimposed on the 0957+561 galaxy counts. Assuming the input galaxy counts 
are correct, the completeness curves are in good agreement, giving us confidence 
that the size- magnitude relationship is close to correct. Fig. shows ellipticity 
histograms for the data and one of the simulations. These are a close match, and 
have similar mean ellipticities of e = 0.265and0.275, respectively. It is important 
that the simulated and real galaxies are well matched in both ellipticity and size 
since the response of a galaxy under shear (the "shear polarizability" - see Kaiser 
et al. 1995) is dependent on these two parameters. The results of the Monte 
Carlo simulations are discussed in the next section. 

6. Discussion 
6.1. 2-d Mass Maps 

The 2-d, KS mass map for our data is shown in Fig. 5 superposed on the 
V-band image of the field. This reconstruction used 1651 galaxies with 24.0 < 
V < 26.5. The detection of the central peak is significant at around the 4.5<r 
level based on the scatter in the simulations. 

The first thing the simulations can reveal are systematic errors in the 2-d mass 
map due to the elongated PSF and the finite size of the image. In Fig. [8] we show 
the average mass map from 102 Monte Carlo simulations of a spherical cluster 
which fits our data well (§ |6.2|) . There is some evidence for small spurious features: 
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Fig. 7. — Ellipticity histograms for galaxies with 24.0 < V < 26.5 in the 0957+561 field 
(solid line) and one of the simulations (dashed line) 24.0 < V < 26.5. The distributions are 
similar, and the mean ellipticities are e = 0.265 and 0.275, respectively. 
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a slight elongation in the north-south direction and la negative features along 
all four edges. These systematics are what is expected from the KS algorithm 
(Schneider 1995), and become more serious as the distance from the image center 
increases. Therefore, in order to quantify measurements made with the mass map 
we will use the simulations as a comparison. 

The isodensity contours in Fig. 5 are not identical to the mean of the simu- 
lations. Is this a significant departure from circular symmetry or just the noise 
in the KS reconstruction? The issue of substructure is an important one from 
the perspective of H determination. There is spectroscopic evidence for the ex- 
istence of a poor cluster at z=0.5 (Garrett et al. 1992, Angonin-Willaime et al. 
1994) centered about 81" west and 43" north of Gl. As pointed out in Bernstein 
et al. (1993), if sufficiently massive, this cluster may be an important constituent 
of the lens. 

In order to test the significance of the deviations from circularity in the KS 
mass reconstruction (Fig. 5), we compute its ellipticity from the quadrupole mo- 
ments of the portions of the map with k > 0. Fig. |] shows this ellipticity value 
along with the distribution of ellipticities measured in an identical manner from 
the Monte-Carlo simulations, which all contain perfectly circular mass distribu- 
tions. The 0957+561 map is a bit rounder than the mean of the simulations but 
well within the normal range. Therefore we conclude that we have not detected 
substructure in our cluster. It should be noted that we have used a smoothing 
scale of s = 30" and substructure finer than this would be washed out. Prob- 
ing higher resolution substructure needs a higher density of background galaxies, 
requiring deeper images with better seeing. 

Also of interest for H Q determination is the accuracy of the mass centroid 
measurement. The mean of the centers of the 102 simulations show no systematic 
deviations from the input center. The la uncertainty in the position of the centers 
is 15". From Fig. 5, the center of mass in the 0957+561 field is located ~ 22" 
from the galaxy Gl (18" east 13" north of Gl). consistent with the position of 
Gl at the 1.5er level. The direction of the offset appears to be consistent with 
that found by Dahle et al. (1994) using 123 galaxies (compared to 1651 in this 
study) although our offset is larger. These two measurements are not completely 
independent since many of the arclets near the QSO will be common to both 
studies. Fig. |10| is a contour plot of galaxy number density for V < 24. There is 
a 3a peak offset from Gl in the same direction as the mass peak lending support 
to the idea that Gl is not at the center of mass. In order to avoid biasing our 
result, we have adopted Gl as the mass center for our azimuthally averaged mass 
density profile described below. 
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Fig. 8. — Mean mass map for 102 simulations of a spherical cluster. The smoothing scale is 
s = 30". The contours have steps of approximately 0.5a of an individual image. Edge effects 
are apparent as well as a small vertical elongation. 
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Fig. 9. — Ellipticy histograms of the simulated mass maps. The input models have e = 
The arrow indicates the measured value for the 0957+561 field. 



Fig. 10. — Contours of galaxy surface number density for V < 24 in steps of 0.5o~ about the 
mean. A total of 302 galaxies were used and a Gaussian smoothing with a = 30" has been 
applied. 
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6.2. Mass Profile 

The azimuthally averaged mass density profile centered on the galaxy Gl as 
derived from the weak lensing approximation (Eqn. [7]) is shown in the upper 
panel of Fig. || The individual points shown in this plot are not independent 
since values at large radii contain a subset of the galaxies used to determine 
quantities at smaller radii. 

Comparing our surface density profile to a variety of spherical cluster simu- 
lations with isothermal plus core profiles we find good agreement between data 
and simulations for: 



S(r) = S [l + (r/r c ) 2 ]- 1/2 , (8) 

with r c = 5" = kpc and So = 3600 ± llOO/i M Q pc~ 2 . There is some 

marginal evidence that the mass profile falls off faster than isothermal for r > 400 
kpc h^ 1 however the uncertainties are large at these radii and a bigger field is 
required to test this hypothesis. 

As can be seen in the upper panel of Fig. |B| the mean surface mass density 
of the simulated clusters determined under the weak lensing approximation is 
slightly steeper than the input surface mass density. This leads to an approxi- 
mately 20% overestimate of the mass at r = 10". This can be attributed to a 
breakdown in the weak lensing approximation. Because we know the true den- 
sity profile for these simulations we can test this by removing the weak lensing 
assumption and seeing if the agreement improves. This amounts to modifying 
Eqn. [7] to the following form: 



7t(r < tv) — 7t(rj < r < r Q ) = — \ - (9) 

K ~ M nt< ro A(iJ)(Ai t 2 + A^) 

— + 

The value of k(R) for each galaxy is calculated assuming the mean S cri t of the 
background galaxy sample. From the Monte-Carlo simulations we find S cr j t = 
5650/i M pc~ 2 . The lower panel of Fig. |B| shows the results of the modified 
mass inversion. The agreement between the input and inferred mass density is 
excellent, with no sign of systematic deviations. The value of the recovery factor 
is < C >= 0.36 ±0.02, very similar to the constant shear HDF result. Therefore, 
for subcritical lenses, one can get accurate density estimates either by fitting the 
data directly to the results of the simulations (which should suffer from the same 
systematics) or, non-parametrically, by adopting an iterative approach where the 
first step assumes the weak lensing approximation and subsequent steps use the 



2 



[l-«(i2)][l-Vl- Di(R) 



-17- 



previous value of the (smoothed) density profile to estimate the jt and (1 — k) 
terms in Eqn. [£| We have adopted the former technique and show the corrected 
profile in the lower panel of Fig. 

Fig. [TT] shows a plot of calibrated < D(r) > in radial bins around Gl. These 
points, though noisy, are independent. The reduced \ 2 between the points and 
model line in Fig. [TT] is 0.62. Due to the large error bars on the points it is not 
possible to provide stong constraints on the shape of the profile. The value of 
< D > for the entire radial range 10" < r < 168" is 0.064 ± 0.019 (s/n = 3.3). 

The Monte-Carlo simulations provide a good check on our error analysis. In 
Fig. ^ we show the surface density and uncertainty derived from the data along 
with the mean and scatter of the simulations. The two sets of error estimates are 
in good agreement. The s/n of the density estimates decreases at the smallest 
radii due to sensitivity to the presence of a small number of radially or tangen- 
tially aligned galaxies. This results because of Ax 2+ Ay 2 in the denominator of 
Eqn. [?]. For example, by removing the four most extreme galaxies, the apparent 
s/n increases from around 3 to 4. 

As mentioned, the best fit model to the data has r c = 5" = llh^ 1 kpc. 
However, this is quite dependent on the innermost point in Fig. || which, as 
just described, is very sensitive to a small number of galaxies. Therefore, it is 
possible that the core radius could be substantially in error. At the low end, it 
is impossible to rule out r c = 0. At the high end, lowering the innermost point 
lcr yields r c = 10" = 33/i _1 kpc. 

One systematic that we do not account for in our Monte-Carlo simulations is 
contamination of our field galaxy sample by cluster galaxies. If cluster galaxies 
have been included in our mass reconstruction then we will have underestimated 
the mass density. If the cluster galaxies are concentrated towards the center of 
the cluster (as one would expect) then we would also be underestimating the 
steepness of the mass profile. Because of the difficulty in distinguishing possible 
cluster galaxies from the field galaxy population we have not attempted to correct 
for contamination. However, a surface number density map of all galaxies used 
in the mass reconstruction (24 < V < 26.5) reveals no significant concentrations 
centered near Gl (or any significant density fluctuations in the field), so we 
conclude that any cluster galaxy contamination is small. 

The surface mass density estimates discussed in this section are determined 
from measurments of density contrasts between innner regions and control an- 
nuli (Eqn. |7|). Therefore an estimate of the total surface density requires that 
we estimate the mass densities of the control annuli. This was implicit in our 
assumption of the functional form in Eqn ^ for the cluster mass distribution. Be- 
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Fig. 11. — A plot of tangential shear vs projected radius in radial bins containing 145 galaxies 
each (10" < r < 168", 24.0 < V < 26.5), centered on Gl. The points are from the data 
and the solid lines are the mean and la upper and lower error bands from 102 Monte-Carlo 
simulations. The dashed line is the input model for the simulations (Eqn. |||) and has r c = 5" 
and Eq = 3600S cr j i . The dashed line is not averaged over the bins. 
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cause this estimate is based on an extrapolation of our measurements to a region 
which we cannot measure, it is worth mentioning the size of the correction. For 
the innermost point (r = 10"), the mean density of the control annulus is around 
7% of the interior mean density if the density profile follows Eqn [5]. This rises to 
about 50% for our outermost point. There is some marginal indication from the 
density profile that the mass density falls off more sharply than our model. If this 
is true then the mass densities of the control annuli are overestimated and the 
mass densities in this paper will be overestimated. Alternatively, if the mass den- 
sities in the control annuli are underestimated (perhaps there is a flat component 
to the cluster mass density), then the mass densities are underestimated. 

7. Other Mass Determinations 

Previous attempts to measure the cluster mass in the 0957+561 field include 
a radial velocity study and an X-ray study. 

Garrett et al (1992) and Angonin-Willaime et al. (1994) measured redshifts 
for 21 probable cluster members, obtaining a velocity dispersion of o = 715 ± 130 
km s _1 . If one assumes a singular isothermal mass distribution and an isotropic 
phase space distribution function, the implied mass contained within a projected 
radius of 1 Mpc is M(r < 1 Mpc) = 3.7±H x 10 14 M . Of course uncertainties 
in the distribution function and problems with field galaxy contamination imply 
a much higher uncertainty. In any case, extrapolating our model with So = 
3600 ± llOO/i M Q pc -2 and r c = YJh~ l kpc implies a mass of M(r < 1 Mpc) 
= 3.9 ± 1.2 x 10 14 M , consistent with the dynamical mass estimate. 

Based on ROSAT PSPC and HRI data of the 0957+561 field, Chartas et al. 
(1995) place upper limits on the cluster temperature and mass. They obtained 
a 3cr upper limit on the cluster mass within a projected radius of 1 Mpc of Gl 
of 1.5 x 10 14 M Q , 2<t lower than the lensing mass estimate. 

8. Conclusion 

The double QSO Q0957+561 was the first multiply gravitationally lensed 
system discovered (Walsh et al. 1979). It is currently the only lensed system 
with a firm time-delay. Extensive modeling has been carried out for this system 
and the conclusion is that in order to constrain Ho the mass distribution in the 
lens must be better understood (Grogin & Narayan 1996 and references therein). 
The lens consists of the galaxy Gl and the cluster to which Gl belongs. 

In this paper we have determined the mass distribution in the z=0.355 cluster 
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in the field 0957+561 by studying the distorted appearance of faint background 
galaxies with 24.0 < V < 26.5. We derived the 2-d mass distribution which 
appears to be consistent with a spherical cluster given our smoothing scale of 
s = 30". The mass centroid is located 22" from the dominant cluster galaxy Gl, 
but is consistent with the galaxy position at the 1.5<r level. A density map of 
galaxies with V < 24.0 has a peak, offset from Gl in the same direction as the 
mass map strengthening the argument that Gl is offset from the center of mass. 

We have calculated the azimuthally averaged mass profile centered on Gl 
out to 400/2T 1 kpc. It is well fit by an isothermal model possessing a core of 
r c = 5" = Uh- 1 kpc. 

In terms of total mass, good agreement is seen with previous kinematic mea- 
surements under the assumption of isotropy. Our gravitational lensing mass, 
however, is substantially higher than upper limits based on X-ray data. 

We will discuss the implications for H in a companion paper. 
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